Brain Topography

Flow, as defined by Mihalyi Csikszentmihalyi (1975), is a holistic sensation experienced when individuals are fully immersed in an activity, resulting in a mental state characterized by a diminished sense of self and altered perception of time. To investigate the global neural dynamics underlying flow, we employed EEG microstate analysis to capture the spatial and temporal properties of dominant transient global brain states (Lehmann et al., 1998). In a study involving 43 participants playing the video game Thumper for 25 minutes, we extracted three four-minute EEG segments from each session corresponding to reported experiences of flow, boredom, and frustration, as determined by self-reports and performance metrics. Across conditions, six distinct microstate topographies (A-F) accounted for most of the global variance. Given that reduced self-referential processing is a key feature of flow, we hypothesized that flow would modulate the properties of microstates C and E, which have been associated with brain regions resembling the default mode network (DMN). Compared to boredom and frustration, the flow condition showed significantly decreased global explained variance, mean duration, time coverage, and occurrence frequency of microstate E, as well as reduced mean duration and time coverage of microstate C. These findings suggest that microstates associated with self-referential processing are shorter and less frequent during flow than during boredom and frustration. This supports the notion that the flow experience modulates global brain dynamics, particularly within the DMN. Furthermore, our results align with previous research reporting reduced DMN activity during meditative and psychedelic states, reinforcing the idea of diminished self-awareness in such conditions.

Sensorimotor beta activity (13-30 Hz) is a key neuronal signature in the human sensorimotor system, and its features can be effectively measured using functional brain imaging methods such as magnetoencephalography (MEG). In addition to its importance in healthy brain processing, beta activity has been shown to be altered in several neurological diseases, underscoring its potential as a biomarker. To serve as biomarkers, features must be reliably defined, stable across measurements and, ideally, amenable to automated analysis, yet current approaches to beta characterization require subjective decisions and manual work. We here describe a hidden Markov model (HMM) based approach to automatically segment beta events from source level MEG beta band activity into discrete high- and low-beta states. We demonstrate the differences between the proposed HMM based approach and a commonly used amplitude-envelope based approach to analyse high- and low-beta modulation. We show that the methods complement each other both when applied to resting data and task related passive movement data. Furthermore, we assess the test-retest reliability of the proposed pipeline within individuals using intraclass correlation coefficients (ICC), and test if HMM constructed at one measurement site can be applied to data acquired at another site, thereby evaluating its multisite transferability. We show that the proposed approach produces stable results within subjects and across sites for many of the features. The ICC values were excellent for high-beta state (86-100% of brain areas), while low-beta state test-retest reliability was more modest. Most of the features showed statistically significant differences between sites only in a few brain areas, indicating very good multisite stability. The proposed approach can serve as an automated, reproducible analysis pipeline for, e.g., clinical applications, and appears suitable for multi-site datasets.

Sex-related differences in the aging of the human brain were studied using large array of experimental data. The open archive CamCan was used as a source of data: the magnetic encephalograms, co-registered with magnetic resonance images of the head, were obtained for each of 434 subjects (ages 18-87 years, mean age 54.7 {+/-}18.4): 217 females (ages 18-87 years, mean age 54.5 {+/-}18.4) and 217 males (ages 18-84 years, mean age 54.8 {+/-}18.3). Recordings were split in 10-year age cohorts, each cohort consisted of equal number of men and women to calculate average intersex characteristics correctly. By massively solving the inverse problem, functional tomograms were calculated - the spatial distribution of elementary spectral components. Physiological noise was eliminated by joint analysis of MEG-based functional tomogram and magnetic resonance image for each subject. Then multichannel spectra were transformed into time series of the power of elementary current dipoles. Summary electric powers were calculated in six conventional frequency bands (1-4 Hz - delta; 4-8 Hz - theta; 8-13 Hz - alpha; 13-21 Hz - beta1; 21-30 Hz - beta2; 30-48 Hz - gamma), and sex differences in age-related changes were examined. It was found that in the youngest age cohort (18-29 years) the summary electrical power of the brain for males is 1.5 times greater than such power for females. For adults (30-69 years), male and female powers are approximately equal, while in older cohorts (70-87 years), male total brain power is greater. Age dependencies in various frequency bands are generally different for men and women, excluding higher frequencies 21-48 Hz. Basic conclusion can be made that after intersex averaging total electric power of the human brain is invariant through the lifespan from 18 to 87 years. The proposed method of joint MEG and MRI analysis can be used for further study of the sex-related details of brain sources in their connection with age changes.

Neural oscillations are fundamental to brain function and cognition. Conventional analyses often rely on predefined frequency bands to assess power modulations, which may obscure finer-grained spectral variability. In this study, we focused on frequency rather than power to investigate whether the natural frequency of each brain region, typically observed at rest, represents a stable intrinsic property or dynamically reconfigures during cognitive processing. We analysed magnetoencephalography (MEG) data from the Human Connectome Project (HCP) across motor execution, working memory, and language processing tasks. Using a multivariate, data-driven spectral clustering approach, we mapped natural frequencies on a voxel-by-voxel basis without imposing predefined bands or regional boundaries. Results indicated that, while the global spatial organization of natural frequencies remained largely preserved during task engagement, specific cortical regions exhibited systematic, task-dependent shifts. In the sensorimotor cortices, the typical resting frequency of [~]24 Hz decreased to [~]6 Hz during movement preparation and at movement onset, and shifted to high-beta rhythms ([~]30 Hz) following hand movement. Increased working memory demands accelerated parieto-occipital alpha/beta activity (from [~]11/16 Hz to [~]13/20 Hz) and recruited high-gamma oscillations (60 to 80 Hz) in medial temporal regions. Finally, arithmetic processing elicited a [~]5 to 15 Hz increase within the beta/gamma ranges across frontoparietal networks relative to semantic comprehension. Taken together, these findings demonstrate that natural frequencies reflect a hybrid architecture: globally stable, yet locally flexible in response to cognitive demands. Moreover, our results suggest that cognitive engagement tends to accelerate neural rhythms in functionally specialized regions, providing a more nuanced understanding of the spectral architecture of human brain function beyond conventional power- and band-based metrics.

Focusing on an organisms task at hand is instrumental for intelligent and goal-driven behavior. However, humans and other animals often fail to pay sustained attention across long time intervals. Failing to stay on-task may cause one to miss crucial task-relevant signals, leading to impaired performance, which can have serious consequences. Therefore, it is important to understand the neural basis of attentional lapses. One promising neural marker of attentional lapses is the frontal P3 (fP3) EEG component, which has been suggested to reflect the susceptibility to incoming sensory input. Following this, we hypothesized that the fP3 1) predicts imminent lapses of attention, and 2) that it should predict upcoming lapses of attention across modalities. In two experiments, we found that the fP3 reliably tracked lapses of attention of sustained attention already seconds preceding the crucial visual signal. We further extended this to the auditory domain: Already 1.5s ahead of the incoming auditory target, the fP3 revealed whether that target was detected or not. Detailed topographic analyses did, however, reveal a slight dissociation between modalities in underlying intracranial source configurations. In sum, this work revealed a supramodal neural signature of susceptibility, which tracks lapses of sustained attention seconds ahead of the critical incoming sensory input.

BackgroundThis study examines a competition-based model (C-model) designed to capture the temporal dynamics of successive brain microstates derived from electroencephalography (EEG) recordings during eyes-open conditions. The analyzed data were obtained from a public repository comprising microstate sequences from 60 sessions of a single subject [1]. When applied to microstate dynamics, the C-model posits a stochastic competition among neural circuits underlying the expression of individual microstates. MethodsThe model is formulated at a conceptual level (computational level in Marrs framework) and employs a geometric distribution to account for the long right tail of microstate duration distributions, interpreted as the probability of "failure" of the currently active microstate to persist. To account for the short-lived left tail, the model incorporates a transient increase in the stability of the currently active network, or equivalently, a temporary decrease in the activation probability of competing microstates (refractory period). ResultsThe model provides a good fit to the microstate duration distributions across all 60 sessions. One third of sessions showed microstate identity sequential dependency with respect to the previous microstates. DiscussionThese results suggest that the C-model captures key aspects of microstate temporal structure. Moreover, because microstate probabilities can be modulated by psychophysiological conditions--including the influence of previously active networks--the model may serve as a building block for more comprehensive neurobiological frameworks of neural and behavioral dynamics. In such frameworks, microstate sequences could emerge from structured competition and flow among neural networks supporting microstate expression.

Conversations can feel effortlessly engaging or, conversely, difficult and unrewarding. Multiple factors contribute to the experienced quality and outcomes of a conversation, among them how interlocutors align with each other. The present study investigated speech-to-speech, brain-to-speech, and brain-to-brain coordination as markers of interpersonal alignment, examining their relationship with jointly perceived interaction quality and mutual affinity between conversational partners. Pairs of previously unacquainted participants (dyads) engaged in multiple short, free-form conversations on topics of varying interest while their vocal and neural activity were simultaneously recorded in a dual-EEG ("hyperscanning") setup. We analyzed interlocutors prosodic adaptation, neural speech tracking, and neural coordination during each conversation. At the speech-to-speech level, our findings reveal that partners with more positive mutual impressions became more similar in their volume and voice quality over the course of the experiment session, reflecting greater prosodic convergence. At the brain-to-speech level, we found no reliable effect of interaction quality on neural tracking of unfolding speech within any individual region, although topographical differences suggested relative modulation across scalp sites. Finally, at the brain-to-brain level, our findings show that higher perceived interaction quality enhanced inter-brain relationships across frequency bands (alpha and theta) and temporal dependencies (concurrent/near-instantaneous and recurrent/listener-lagging), with the strongest effects observed for concurrent alpha-band coupling. These findings suggest that distinct coordination processes are involved in how interlocutors experience an interaction and how they establish relational affinity, casting new light into the mechanisms that make a conversation worthwhile.

The routing of information across the brains structural network is central to its wide range of functional capabilities. However, the mechanisms underlying information routing in complex brain networks, particularly between regions that do not share a direct anatomical connection, remain poorly understood. Neural mass models (NMMs), a computational modelling framework capable of capturing complex neural dynamics across scales, can potentially be used to study the dynamical and network bases of these vital polysynaptic routing processes. In this study, we investigate polysynaptic signalling in three widely used NMMs, obeying Ornstein-Uhlenbeck, Stuart-Landau, and Jansen-Rit dynamics, by tracking the propagation of a discrete, focal, high-amplitude perturbation across the underlying network. We find that polysynaptic propagation emerges in all tested NMMs when configured within dynamical regimes that effectively enhance the persistence of perturbations. We also find distinct parameter domains that maximise signal propagation to directly connected regions and to those separated from the source by at least two hops. Finally, we benchmark in silico stimulus propagation in the brain network against an empirical dataset of direct electrical stimulation trials, to explore the relative capabilities of the NMMs in capturing signal propagation to connected versus unconnected regions. This analysis highlights the significance of dynamical repertoire in capturing stimulus propagation outcomes. Overall, this study provides insights into how dynamical and network features shape signal propagation over complex brain networks.

Vocal communication involves a series of cognitive processes, which can be broadly categorized into three components: perceiving communicative signals; deciding whether and how to respond; and generating vocal motor output. These processes must work harmoniously, with integration and bridging between components being crucial for effective communication. Previous research on vocal communication has typically focused on specific brain regions or isolated cognitive functions, often lacking a holistic perspective of macro-scale, whole-cortical dynamics and their role in the complete communication process. Therefore, although the cortical areas associated with each cognitive component have been localized in humans, the macro-scale cortical dynamics underlying the integration of these cognitive processes remain unknown. Building on recent findings linking macro-scale cortical dynamics to behavioral performance, we hypothesized that traveling wave like cross-areal interactions play a role in integrating the three communicative components. To test this hypothesis, we recorded whole-cortical activity using epidural electrocorticography (ECoG) while subject marmosets vocally interacted with partners. We found theta-band activation in several cortical areas, including the parietal and auditory cortices, while listening to partners calls. This activity was further modulated depending on whether the subjects engaged in vocal interactions, potentially representing the transformation of sensory processing into decision-making and vocal motor preparation. Given the widespread nature of this modulation, we next characterized whole-brain activity patterns by employing a novel analytical method, Weakly Orthogonal Conjugate Contrast Analysis (WOCCA). This analysis revealed that cortical activity could be decomposed into two distinct traveling wave like propagation patterns, a rotational and a translational wave, and both waves discriminated communicative conditions consistent with localized activity. The rotational wave further represented vocal motor preparation through trigger-like temporal pattern. In addition, the magnitude of the translational wave immediately before subjects vocal production correlated with the vocal production-induced suppression of high-gamma-band activity, particularly in the prefrontal and auditory cortices. As vocalization-induced suppression is believed to reflect sensory prediction, the translational wave may propagate specific decision-related or acoustic information necessary for subsequent vocal production to local cortical areas. These findings suggest that the brain orchestrates the sequential cognitive processes underlying vocal communication through macro-scale traveling waves.

Most of what we currently know about brain oscillations is derived from Fourier-based spectral power methods. While widely adopted, these procedures introduce methodological limitations and inherently overlook informative neurophysiological features that often remain unreported. In this study, we characterized resting-state oscillatory activity from human magnetoencephalography (MEG) recordings (N = 128) by integrating two complementary approaches. First, oscillatory episodes were detected at the source level with sBOSC. Subsequently, the ByCycle algorithm was applied to these episodes to extract individual cycle features. Results revealed that the brain engages in oscillatory activity for only [~]25% of the recording time, with occipital and parietal regions accounting for the highest temporal prevalence across canonical frequency bands. Furthermore, oscillatory episodes lasted an average of 4.6 cycles, reinforcing the view of neural oscillations as transient bursts. Region-specific duration and power measures revealed distinct anatomical organizations offering complementary physiological information. Finally, by extracting the instantaneous amplitude, period, and waveform asymmetry of individual cycles, we successfully dissociated sinusoidal occipital alpha waves from the asymmetric sensorimotor mu rhythm. By moving beyond traditional power-centric analyses, this approach provides a comprehensive characterization of spontaneous oscillatory activity, thereby offering new insights into the spatial, temporal, and spectral structure of human brain oscillations.

The heartbeat-evoked potential (HEP) reflects the cortical processing of cardiac afferent signals. However, it remains unclear whether trial-level interoceptive prediction errors can be quantified directly from spontaneous resting cardiac fluctuations and whether these model-derived errors are associated with HEP amplitude. Here, we applied a Kalman filter, implemented as a sequential Bayesian estimation procedure, to resting-state EEG and ECG recordings from 21 healthy adults to estimate trial-by-trial signed prediction errors in RR-intervals. Positive prediction errors reflected unexpected cardiac deceleration, whereas negative prediction errors reflected unexpected cardiac acceleration. Cluster-based permutation tests showed that unexpected cardiac acceleration was associated with greater fronto-centro-parietal HEP amplitude than unexpected deceleration in an early post-R-peak window, spanning FC1, CP1, Pz, CP2, Cz, C4 and FC2 from 215 to 250 ms. A Bayesian linear mixed-effects model further indicated a credible negative association between signed prediction error and HEP amplitude after controlling for respiratory phase and preceding RR interval. In a secondary connectivity analysis, unexpected acceleration was associated with stronger Cz-to-frontal beta-band phase synchrony during a later post-R-peak window from 250 to 500 ms. Exploratory individual-difference analyses suggested that neuroticism was negatively correlated with late frontal HEP amplitude during unexpected acceleration, but not during unexpected deceleration or when trials were pooled across conditions. These findings demonstrate that spontaneous cardiac fluctuations can be used to derive trial-level computational estimates of interoceptive prediction error and that these estimates are reflected in early HEP amplitude. They further suggest that the cortical processing of unexpected cardiac acceleration may be related to individual differences in affective personality traits.

The awake brain is known to display asynchronous (AS) states during periods of attention and arousal, but the responsiveness properties of such states remain unclear. Here, we investigate this question using computational models of spiking networks of excitatory and inhibitory neurons, mimicking recurrently-connected networks in layer 2/3 of the cerebral cortex. The networks can generate a continuum of AS states, but with different responsiveness characteristics. By using a mean-field model to infer the dynamic properties of the system, we find that there are two families of AS states, which we call "underdamped" (UD) and "overdamped" (OD). Responsiveness is maximised at the transition between OD and UD states, which identifies a "working point" that may present advantageous computational properties.

Growing interest to describe the electrical behavior of glial cells, mainly astrocytes, in intact brain tissue poses more and more challenges to commonly accepted belief they only respond in a linear manner in uptake of the excess of extracellular potassium and maintenance of their network equipotentiality. Their highly conductive mutual interconnections via gap junction (GJ) connections introduce yet another class of nonlinear elements. As more studies report nonlinearities in membrane voltage Vm dependence of both, the membrane and junctional conductances, the need to formulate minimal dynamical models of their transient behavior is getting more acute. Since ODE models of coupled cells, even in simplest 1-d arrays, require simplified descriptions and small set of parameters, rare quantitative studies on glia makes the task even more difficult. This study attempts to qualify a self-coupled cell, or a glial cell coupled to fixed voltage as useful system for detecting the nature of instabilities and transitions coming from coupling. In a novel biophysical model of coupled astrocyte, we introduce nonlinear kinetics of deactivation for large junctional voltages for the first time. We found that N-shaped nonlinearities and corresponding fold structure in the vector field of isolated cell serves as a baseline on top of which coupling nonlinearities enrich the bifurcation picture. Numerical simulations of 1-d array of coupled astrocytes show that coupling increases the propensity of astrocytic Vm to bistability and front propagation. We believe that presented illustrations of possible effects of coupling nonlinearities will motivate neurobiologists to further explore their impact in disease. Significance statementTransient changes in membrane voltage of glial cells may produce significant transient voltage difference between directly coupled cells. Nonlinear steady-state conductance of their interconnection elements, the gap junctions, introduce nonlinear current profiles which are very difficult to measure and quantitate using the available methods due to marked permeability of the junctions and leakiness of glial membrane in general. We propose a minimal model of glial membrane extended with a self-coupled feedback loop, which under realistic simplifying assumptions could serve for qualitative analysis of the impact of coupling, on the stability of resting membrane voltage. Neuronal cells of the brain and spinal cord cannot exist and function without supportive and neuromodulatory functions of the diverse population of glial cells. This applies to virtually all physiological processes on cell level - from cell development, metabolic support, membrane signaling, slow molecular signal transduction, ion homeostasis, neurovascular coupling, myelination, to mention only a few, manifest neuro-glial interaction. Even though all glial cell types are interconnected, the most abundant ones, the astrocytes are massively interconnected by gap junctions to form ordered networks. Electrically, astrocytic networks display membrane voltage equipotentiality, which is considered system-wide resting state for given neuro-glial circuit or unit. With molecular and cellular substrates of glial connectivity being slowly elucidated, network science and dynamical modeling are slowly "invading" that area with many important issues left open. In this study using classical dynamical systems approaches we give indications how nonlinear intercellular coupling between astrocytes affects physiological resting state and its instabilities compared to isolated, uncoupled cell. We strongly believe the suggested minimal model could fill the gap in ODE modeling of neuro-glial circuits, within broadest scope of hypothesis-driven research in cell-level neuroscience.

Pig vocalizations convey information about the emotional states of individuals, varying with arousal and valence. Studies show that different call types reflect distinct emotional contexts and social interactions for the receivers. However, little is known about the brain mechanisms behind the perception of conspecifics vocalizations. This study used BOLD fMRI to explore how pigs brains respond to emotionally varied vocalizations, with the aim to identify activity in regions linked to emotion, reward, and social processing. Eight healthy 2-month-old pigs underwent auditory brainstem response (ABR) testing and BOLD fMRI to assess brain responses to pig vocalizations with different hedonic valence. Sounds were delivered via MRI-compatible earphones, and imaging was performed on a 1.5T scanner. Data were analyzed using voxel-based and ROI-based statistics in SPM12 with small volume correction (SVC). Due to hearing anomalies or MRI artefacts, only 5 pigs were included in the final analysis. Functional MRI revealed that vocalizations activated regions of the auditory pathway and the left amygdala (pFWE at peak < 0.05 after SVC for all), with specific differences between positive and negative sounds. Clusters of activated voxels covering part of hippocampal areas, caudate nuclei and putamen were found with both positive and aversive vocal sounds. Limbic regions, including the amygdala and insula (p<0.05), as well as the right hippocampus after SVC (pFWE = 0.015) were uniquely engaged during the perception of negative conspecific vocalizations, indicating distinct processing based on emotional valence. This study shows for the first time that piglets brain can process and differentiate emotional vocalizations from other pigs, even under general anesthesia. Positive and negative vocal sound playbacks activated distinct brain regions related to hearing, emotion and reward. These findings highlight pigs cognitive and emotional processing of vocal cues. This study is part of a wider research program aimed at developing the fMRI protocol with acoustic stimulation in juvenile pigs.

IntroductionStimulation-evoked potentials (SEPs), recorded both during and after deep brain stimulation (DBS) surgery, have shown promise for guiding DBS targeting and programming. However, filtering protocols applied to stimulation trains produce an artifact we call a filter-induced oscillation (FIO) which closely mimics physiological SEPs. Hence, we outline the mechanistic origins of this distortion and describe a means of differentiating it from valid SEP activity. MethodsWe recorded in 18 patients undergoing DBS surgery targeting the subthalamic nucleus or globus pallidus internus. We stimulated target nuclei with cathode-first (CF) and anode-first (AF) pulses to record native SEPs, and in white matter tracts (null condition). Recordings were subsequently filtered to illustrate FIO. Next, we filtered harmonic frequencies of an artificial stimulation train to demonstrate FIO origins. Finally, FIO was deliberately generated in white matter recordings with a notch filter, and its behaviour contrasted with SEPs during AF and CF stimulation. ResultsFiltering stimulation trains produced FIOs that depended on filter order and corner frequency. We also showed that FIO emerges from filter-induced attenuations of harmonic frequencies which compose stimulation trains, producing oscillations of like frequency around pulses. Finally, FIOs reverse in polarity depending on AF or CF stimulation, whereas SEPs do not. ConclusionsGiven the potential for widespread adoption of SEPs in DBS targeting and programming, safe analytical protocols are imperative to avoid the induction of processing-related artifacts which can be misinterpreted as biological signals. Here we establish the necessary theory for identifying FIOs and tuning analytical pipelines to avoid their generation.

Transcranial temporal interference stimulation (tTIS) is a promising non-invasive brain stimulation technique that has the potential to selectively modulate deep brain regions by delivering two high-frequency alternating currents that interfere to produce a low-frequency amplitude-modulated envelope at the target. A key challenge in deploying tTIS is finding electrode current patterns that are simultaneously effective, focal, and safe. This is a fundamentally non-convex optimization problem for which a number of methods have recently been proposed. However, no systematic comparison of these methods across a large and diverse set of brain targets has been performed, leaving practitioners without clear guidance on how best to optimize for a particular experimental setting. In this work, we present a comprehensive benchmarking study comparing seven tTIS optimization methods that have appeared in the literature in recent years: exhaustive search, genetic algorithm, multi-objective evolutionary algorithm (MOVEA), min-max optimization, convex TI (CVXTI), non-convex optimization with convex relaxations, and an unsupervised neural network. All methods were evaluated across 250 brain targets spanning cortical and subcortical gray matter and white matter regions in five finite element head models. Each method was evaluated on two key metrics: mean electric field strength within the target region of interest, and off-target stimulated brain volume. Results were stratified by tissue type and target depth to identify systematic performance differences. Based on these results, we provide practical evidence-based recommendations for optimization method selection among these seven methods depending on target location, tissue type, and available computation time. Moreover we provide the code base that will allow other investigators to use these methods for their own applications. Our goal is to provide researchers and clinicians with a clear, evidence-based framework for choosing a tTIS optimization method suited to their specific target and application.

Pleasurable urge to move to music is often referred to as groove. Although previous studies have shown an association between the supplementary motor area (SMA) and the groove experience, its causal role remains unclear. Here, we investigated whether the SMA is causally involved in groove experience during music listening using repetitive transcranial magnetic stimulation. Fifteen healthy participants completed three sessions on separate days: excitatory stimulation (intermittent theta burst stimulation; iTBS) over the SMA, inhibitory stimulation (continuous theta burst stimulation; cTBS) over the SMA, and sham stimulation (iTBS or cTBS) over the vertex. After each stimulation session, participants listened to five high-groove and five low-groove musical excerpts and rated urge-to-move and pleasure on a 0-100 scale. Heart rate was additionally recorded as an exploratory physiological measure during music listening. Linear mixed-effects models (LMM) showed that pleasure ratings, but not urge-to-move ratings, were higher following both iTBS and cTBS compared with sham stimulation. In exploratory LMMs, reduced log-transformed heart rate variability (HRV) significantly predicted higher pleasure ratings. These findings suggest that SMA stimulation modulates the pleasurable component of the groove experience, likely via network-level mechanisms rather than a simple linear relationship between SMA excitability and pleasure. They also raise the possibility that reduced parasympathetic activity, reflected by lower HRV, mediates the stimulation-related increase in musical pleasure. Future studies should investigate the causal roles of other brain regions as well as clarify the directionality between autonomic changes and the groove experience.

Functional connectivity in network neuroscience is traditionally characterized using time-averaged correlations between brain regions. While these summaries capture stable large-scale organization, they do not fully reflect the temporal structure of moment-to-moment interactions. Here, we investigate how the order of interaction used to represent brain dynamics shapes the organization recovered from neural data. We compare three interaction representations of fMRI dynamics: regional activation (node time series), pairwise co-fluctuations (edge time series), and higher-order triplet interactions (triangle time series); within a common topological framework using Mapper from topological data analysis (TDA). Across task and resting-state data, Mapper representations derived from pairwise co-fluctuations more distinctly segregate task conditions than activation-based or higher-order representations. This organization reflects structured coordination patterns beyond activation polarity and is driven by high-amplitude interaction events. Beyond task states, modularity quality computed across all Mapper representations is highest for edge time series and selectively associated with stable individual differences: higher modularity relates to higher conscientiousness and lower internalizing and externalizing symptom dimensions. Together, these findings suggest that behaviorally relevant information is reflected in the topology of moment-to-moment brain interactions. Topological analysis of interaction-level dynamics therefore provides a complementary and interpretable framework for linking large-scale neural coordination to cognition, personality, and mental health.

Objective. Continuous, unsupervised monitoring of cognitive brain responses has long been constrained by the demands of laboratory EEG. Whether the P300 event-related potential, an established marker of attention and cognitive processing, can be elicited as an incidental byproduct of genuine gameplay, recorded with a minimal wearable EEG system under unsupervised home conditions, has not been established. Approach. Ten healthy adults played a gamified visual oddball task in which infrequent target stimuli (green gates) were embedded among frequent non-targets (red gates) within a continuous third-person running game. EEG was recorded with a four-channel dry-electrode headband (EEG channels: O1, O2, T3, T4; forehead reference; 250Hz) with self-mounted electrodes in a home setting, without experimenter supervision. Group-level effects were assessed with cluster-based permutation tests and peak-amplitude tests. Single-trial classification used linear discriminant analysis (LDA) with four features per channel (16 total). Additional analyses included a within-subject comparison with a classical visual oddball paradigm using identical hardware, pilot data from a patient with relapsing-remitting multiple sclerosis, within-subject stability across 48 sessions, and pilot recordings with a headphone form factor. Main results. A robust P300-like difference wave emerged on all four channels at the group level (cluster-based permutation tests, p < 0.05), with individual-level detection in 8 of 10 participants (exact binomial p < 0.001). Single-trial LDA yielded a median cross-validated AUC of 0.730 (95% CI 0.672-0.820), with 9 of 10 participants exceeding chance. In a within-subject comparison, waveform morphology was closely preserved relative to a classical laboratory oddball, and classification performance was markedly higher in the game condition (AUC 0.820 versus 0.555). A patient with relapsing-remitting multiple sclerosis produced a clear P300 (AUC 0.853) with latencies within the healthy range. Within-session split-half reliability was high (r > 0.70 on three of four channels), though between-session reliability was near zero across 48 sessions in one participant, with a declining classification trend over time. Pilot recordings with a headphone form factor also yielded a P300-like deflection. Significance. These results demonstrate that the P300 can be elicited as a gameplay-integrated neural readout during genuine gameplay with a wearable, dry-electrode EEG system under unsupervised conditions. Gamification does not compromise P300 elicitation; in the within-subject comparison, it enhanced single-trial discriminability. The findings indicate that gamified, home-based P300 monitoring is achievable with minimal hardware and provide preliminary evidence for applicability in clinical populations, most notably multiple sclerosis, where P300 has established biomarker value but where the logistical burden of laboratory assessment currently precludes longitudinal use.

The sensory content and temporal structure of stimuli have been shown to consistently bias duration perception. Temporal intervals filled with continuous sensory input ("filled intervals"), are often perceived as lasting longer than intervals marked only by their onset and offset ("empty intervals"). Despite this robust behavioral finding, it remains unclear whether filled and empty intervals rely on similar or distinct neural mechanisms and, more generally, how sensory format shapes the neural processing of millisecond time. To address this question, we asked twenty-one healthy participants to reproduce visual durations across different stimulus configurations while high-density scalp EEG was recorded. Behavioral results revealed differences in performance across stimulus configurations. Event-related potentials (ERPs) recorded at occipito-parietal and fronto-central electrodes between 0.1 and 0.4 s after duration offset were modulated in amplitude by both stimulus duration and format. These modulations scaled with the sensory load of the stimulus and its duration, suggesting a common underlying mechanism. A Representational Similarity Analysis (RSA) of the ERP data showed that perceived time was represented more strongly than physical time particularly at occipito-parietal electrodes, but only within the 0.2-0.3 s post-offset window, where stimulus format exerted a pronounced effect on the ERP signal. These findings highlight the role of sensory processing in shaping duration perception and its neural coding, and reveal an early neural signature of perceived time in occipito-parietal electrodes. 1 Significance statementOur perception of subsecond durations is distorted by the sensory content of stimuli. Here, we investigated how stimulus configuration shapes the neural correlates of visual duration perception. Specifically, we asked whether temporal intervals filled with continuous sensory input are processed differently from those lacking such content. We found that, between 0.2 and 0.3 s after interval offset, ERP amplitudes were modulated by stimulus content, and in this same temporal window the EEG signal reflected the perceptual bias. These findings support the view that duration processing and perception are deeply rooted in sensory processing.